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ABSTRACT 


Heat exchangers with one shell pass and n tube passes 
are often referred to as 1l-n exchangers. The heat transfer 
literature contains many references to studies of l1-n 
exchangers when n is even but apparently, other than a 
single study and some work pertaining to the 1-3 exchanger 
at the Naval Postgraduate School, little work has been done 
with respect to the l-n exchanger when nis odd. This 
thesis expands upon the study of l1-n exchangers with n being 
odd. 

While a completely closed form solution was found to be 
unfeasible, a polynomial approximation has been developed 
that yields the effectiveness (ce) of the two possible 
arrangements of the 1-5 exchanger as a function of the 
capacity rate ratio (R) and the number of transfer units 
(N+,)- This will enable the analyst to consider exchangers 
where the inlet to and outlet from the tubes are at opposite 


ends of the exchanger. 
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NOMENCLATURE 


English Letter Symbols 


A4 


As 


Exchanger heat-transfer surface, sq m 
Coefficient of the mth value dimensionless 
Coefficient, dimensionless 


Ist order coefficient to be multiplied by Ney > 
dimensionless 


2nd order coefficient to be multiplied by Nyy, 
dimensionless 


3rd order coefficient to be multiplied by Ny,3, 
dimensionless 


4th order coefficient to be multiplied by Neu; 
dimensionless 


5th order coefficient to be multiplied by Ny,,°, 
dimensionless 


Exchanger heat-transfer surface, sq m/m 


Capacity rate, W/°K. Also designates dimensionless 
arbitrary constant 


Specific heat at constant pressure of cold fluid, 
J/kg°K 


Specific heat at constant pressure of hot fluid, 
epics — 1K 


Empirical value of effectiveness (computer generated), 
dimensionless 


Error. Also used as the exponential function 


Logarithmic mean temperature difference correction 
factor, dimensionless 


Exchanger length, m 


Number of effectiveness empirical data points used to 
determine a curve for R, dimensionless 


Number of tube passes, 
equations, 


dimensionless. 


dimensionless 


Also number of 


A related Nyy per unit length cold side, m71 


A related Nt, per unit length hot side, m7! 


Number of transfer units, 


Temperature group, 


dimensionless 


Total rate of heat transfer, W 


Capacity rate ratio, 
Temperature group, 


Sum of the squares of the residuals, 


= Maximum total rate of heat transfer, 


Hot fluid temperature, 


Particular integral, 
Hot fluid temperature in, 


Hot fluid temperature out, 


Cold 
Cold 
Cold 
Cold 
Cold 
Cold 
Cold 
Cold 
Cold 
Cold 


Cold 


‘Ee ILiol aks! 
fluid 
fluid 
fluid 
Elude 
fluid 
£ Tone 
fluid 
EA MEL alco! 
£ Liblauel 


fluid 


temperature 
temperature 
temperature 
temperature 
temperature 
temperature 
temperature 
temperature 
temperature 
temperature 


temperature 


dimensionless 


dimensionless 


rc 


dimensionless 


ake 
ae 
sly oh rutile © 
Out. © 
lst pass, °C 
between ist and 
2nd pass, °C 
between 2nd and 
Sra pass suse 
between 3rd and 
4th pass, °C 
between 4th and 


Sth pass, °C 
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dimensionless 


W 


dimensionless 


2nd passes 


3rd passes 


4th passes 


oth passes 


U = Overall heat transfer coefficient, W/m? - °C 

W = Mass flow, kg/sec. 

x = length coordinate, m. Also used to represent 
a constant value in a sequence 

y = Sum of mth degree polynomial, defined by eq. (5.1), 
dimensionless 

Greek Letter Symbols 

at = Root of auxiliary differential equation, 1/m 

€ = Exchanger effectiveness, dimensionless 

x = Combination of variables defined by equation (3.12), 
dimensionless 

» = Summation, dimensionless 

o2 = Variance, dimensionless 

6, = Mean temperature difference for exchanger, °C 

) = Indicates partial derivative, dimensionless 

8 = A combination of terms defined by eq. (3.23) 

Subscripts 

C = Cold fluid 

h = Hot fluid 

i,j,k = Values in a sequence 

m = Degree or order, an exponent 

1 = inlet 

2 = outlet 
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I. INTRODUCTION 


In many engineering applications, it is desirable to 
transfer thermal energy between fluids at different 
temperatures. This action may be economically accomplished 
with the aid of a heat exchanger. 

The primary function of a non-direct contact heat 
exchanger is to provide two paths of flow, one for the hot 
fluid and one for the cold, by means of which heat can be 
transferred through walls separating the fluids. 

The quantity of heat transferred is governed by three 
factors: 1) the amount and nature of the heat transfer 
surface to which the two fluids are exposed, 2) the overall 
heat transfer coefficient, and 3) the mean temperature 
difference across the intervening wall from one fluid to the 
other. 

Numerous flow arrangements and geometric designs are 
possible. One of the most common types of heat exchangers 
is the shell and tube type. Here one fluid flows through 
the tubes. The tubes are enclosed in a Shell with 
provisions for the other fluid to flow through the spaces 
between and exterior to the tubes. Overall flow may be 
either parallel or counter, but in most cases, the latter is 
preferred. The reason for this will be seen later. One of 


the most common flow arrangements used in the Shell and tube 
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heat exchangers is a combination of the two, the parallel=- 
counterflow exchanger. Heat exchangers with one shell pass 
and n tube passes are therefore referred to as l1-n heat 
exchangers. 

One example of this flow arrangement is the one 


shell-pass two tube pass exchanger as shown in Figure Taam 


+ 


1 
\ Baffle 
Qt =< 


= 


=~ 


Figure 1.1 A 1-2 Heat Exchanger (One Shell Pass and Two 
Tube Passes) 
This type of exchanger is configured so that all of the tube 
side fluid flows consecutively through half of the total 
number of tubes: The baffles in the shell serve to sib l>-< 
the shell fluid and to maintain fluid flow at right angles 
to the tube passes in order to obtain a maximum "shell 


side" heat transfer coefficient. 
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O'Hare [Ref. 1] has indicated that, to date, much work 
has been done on finding the logarithmic mean temperature 
difference for heat exchangers with an even number of passes 
but little has been done with exchangers having an odd 
number of passes. From this basic statement, he proceeded 
with a theoretical examination of the 1-3 heat exchanger 
which resulted in some startling conclusions. As a follow 
up to his work, this thesis will parallel and follow the 
format of O'Hare [{Ref. 1] for the 1-5 parallel counter flow 
exchanger. One arrangement of the 1-5 exchanger is 


displayed in Figure 2. 


- 
{7 Baffle 
a a eh 

| a aa 


Figure 1.2 A 1-5 Heat Exchanger (One Shell Pass and Five 
Tube Passes) 


ats, 


II. THE DEVELOPMENT OF THE EFFECTIVENESS METHOD 


A. HISTORICAL DEVELOPMENT 

Nagle [Ref. 2: pp. 604-609], in 1931, credited Davis 
[Ref. 3] with a simplified method for computing actual temp- 
erature differences between two heat-interchanging streams 
which depart from true counter or concurrent (parallel) 
flow. This is now the familiar "F factor" method which 
expresses the actual mean temperature difference 8, in q = 
UA8, as a fraction F of the counterflow logarithmic mean 
temperature difference, LTIMD, 9,6 via 9m = FOme- 

The example of initial interest was the 1-2 exchanger 
with a single shell pass and two continuous tube passes in 
counter and concurrent flow with it. The method involved 
derivation of the actual temperature difference for the flow 
pattern and formed the ratio F = 07)/®mc- This familiar LMTD 
correction factor was plotted conveniently as functions of 
the effectiveness, ¢«, and the capacity rate ratio R with R 
aS a parameter. These mean temperature difference 
correction charts are available for many flow arrangements 
[Ref. 4: pp. 829-833 and Ref. 5]. The effectiveness, e 
(often called P or S), is always the cold fluid 
effectiveness and R is always the capacity rate ratio of 


eold fluid £o Hot ft lhiume- 
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Nagle detailed assumptions and derivations for the 1-2, 
1-4 and 1-6 exchangers. The F factors were obtained by 
Nagle through graphical integration and were accompanied by 
the comment that F factors for the 1-2 exchanger could be 
applied with negligible error to 1-4 and 1-6 exchangers. 
Underwood [Ref. 6: pp. 145-148] rederived the equations of 
Nagle for 1-2 and 1-4 exchangers to eliminate the need for 
Obtaining F factors by graphical integration. 

Bowman [Ref. 7: pp. 541-544] pointed out that for a 
very large or infinite number of tube passes, the F factor 
approached, as a limit, its value in crossflow with both 
fluids completely mixed. [It was further stated that even at 
the limit, the F factors were only 1 to 2 percent lower than 
those for the 1-2 exchanger. A previous paper by Kraus and 
Kern [{Ref. 8] did not confirm the generalization that l-n 
exchangers differed only negligibly from the 1-2 exchanger 
although this lack of confirmation was obtained on an 
e = £(R, Ny) basis. Moreover, the Kraus-Kern work does not 
confirm the generalizations on an F = f(R, Ney, €) basis. 

From the standpoint of usefulness and good accuracy, it 
is essential that F factors, if they are to be used in 
preference to e« = f(R, Ney, flow arrangement), be obtained 
with precision. Plots of F = f{(R, e«, = P or S) [Ref. 4: 
pp. 829-833 and Ref. 5] show that the curves for particular 
values of R approach infinite slope as F decreases. While 


this can be partially alleviated by restricting R < 1.0 (a 


ir 


constraint used in the « = {(R, Ny,, flow arrangement) 
approach), it is seen that small errors in the interpolation 
for R or e = P or S can result in large fluctuations in the 
value of F. 

In a comprehensive paper, Bowman, Mueller and Nagle 
[Ref. 9: pp. 283-294] presented graphs of F factors for 
shells with one through six shell passes and numbers of 
continuous tube passes respectively double the number of 
shell passes. In view of the earlier references to Nagle 
and Bowman, it should be noted that F factors were computed 
for the 1-2 exchanger in [Ref. 9: pp. 283-294] using the 
equations of Underwood [Ref. 6: pp. 145-158]. 

Ten Broeck [Ref. 10: pp. 1041-1042] prepared a graph 
of the dimensionless groups now known as ec, R and Ny, for 
the 1-2 exchanger. Such a graph had the added versatility 
of simplifying the calculation of performance in a given 
exchanger when operating at conditions different for those 
for which it was designed. Kays and London [Ref. 11: 
pp. 63-74] prepared similar graphs and tables of 
e = f(R, Ney, flow arrangement) for the 1-2 exchanger and 
for several cases of crossflow and periodic flow. 

In the foregoing discussion, it has been observed that 
there are two methods for the design and analysis of a heat 


exchanger. These are the so-called F factor” and eai7a 
approaches and it is to be noted that the Heat Transfer 


Literature contains many references of design equations for 
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either of these methods for exchangers that have one shell 
pass and even numbers of tube passes. 

It was O'Hare [Ref. 1] who pointed out that there may be 
a virtue in having an odd number of tube passes particularly 
in a marine application. Such an arrangement allows the 
tube side fluid to enter one side of the exchanger and leave 
from the opposite side (see Figure 1.2). 

With the possible use for an odd tube pass heat 
exchanger established, a literature Search was conducted to 
uncover any previously solved cases of l-n heat exchangers 
with n, the number of odd tube passes. O'Hare [Ref. 1], 
Fischer [Ref. 12], and Stevens, et. al. [Ref. 13] have 
dealt with this type of exchanger. Although Stevens, 
et. al. deal with multiple crossflow passes, all provide the 
foundation for the work reported here as well as additional 
work pertaining to the 1-3 exchanger. O'Hare's investiga- 
tion centered around the development of a solution for the 
Beoeexchanger in the effectiveness-N;,, framework while 
Fischer derived an equation using the mean temperature 
difference for the 1-3 exchanger. Moreover, Fischer 
expressed his results in terms of the correction factor F 
and treats only the case of the 1-3 exchanger. This thesis 
will expand on the N;,, method of thermal heat exchanger 
design by covering the 1-5 heat exchanger. The work will 


treat both the three counterflow and two parallel flow 
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and the two counterflow and three parallel flow (1-5: 3P) 


cases (see Figures 2.1 and 2.2). 
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Figure 2.1 1-5:3C Five Tube Passes - Three in Counterflow 


Breure 2.2 1-5:3P Five Tube Passes - Three in Parallel Flow 
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III. AN ATTEMPT AT AN ANALYTIC SOLUTION 
A. EFFECTIVENESS AS A FUNCTION OF CAPACITY RATE RATIO AND 

EXCHANGER SIZE (Nty) 

This section deals with an investigation into the 
effectiveness, e«, of a one shell pass-five tube pass heat 
exchanger, in which e« compares the actual heat transfer 
capability to the thermodynamically limited, maximum possible 
heat transfer capability. This exchanger heat transfer 


effectiveness is given by 


CaCr OG 


hot,in Thot,out? me ec 


) 


cold,out | vGold,in 


E 


max min‘ hot,in cold,in min ‘ hot,in co ldyarna® 


where C,;, is the smaller of the Cy, and C, magnitudes (the 
capacity rates). Thus, ¢ possesses the significance of 
effectiveness of the heat exchanger from a thermodynamic 
point of view, with the magnitude of the effectiveness 
completely defining the heat transfer performance. In 
general we express e« = f(Nty, R, and flow arrangement) and 
when the flow arrangement is specified or understood, it is 
said that « = f(Ney, BR). [Ref. 11: pp. 14-26]. 

The number of heat transfer units Nt, is a nondimensional 
expression of the "heat transfer size" of the exchanger. 
When N+, is small the exchanger effectiveness is low, and 


when N+, is large, « asymptotically approaches the limit 
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imposed by the flow arrangement and thermodynamic 


conditions. From inspection of the definition of My 


A 
ee eee ol ‘ 
Mut ee rc f ete 
min min % 
it is clear that the overall conductance and transfer 


area affect the costs of attaining a high value for Ny, 


ergo a high «. The capacity rate ratio, R, as defined by 


is simply the ratio of mass flow rate times specific heat 
capacity for the two streams. These can be considered as 
flow stream thermal-capacity rates, l.e., tee energy storage 
rate in the stream per unit of temperature change. [Ref. 
ime pp. 14-26] 

The attempt taken in this thesis to develop a closed 
form solution has used the basic fundamentals of heat 
transfer as well as the foregoing definition. A closed form 
solution for ¢« was sought for both 1-5 exchangers with one 
having three out of five tube passes in parallel flow and 
the other having three out of five tube passes in counter- 
flow. The analytical approach taken, and demonstrated in 
this section, is for three out of five tube passes in 


counterflow. 
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B. ANALYTICAL DEVELOPMENT 

The derivation for the effectiveness, ce, of the 1-5 
exchanger as a function of the capacity rate ratio, R, and 
number of transfer units, Ny,, depends on several 
assumptions. 


(1) The overall coefficient of heat transfer, U, does 
not vary within the exchanger. 


(2) The specific heat of both hot side and cold side 
fluids does not vary. 


(3) Each fluid is thoroughly mixed, that is, the 
temperature of both hot and.cold side fluids is 
uniform over any cross section. 


(4) Steady flow conditions are maintained. 


(5) Heat losses to or from the environment are 
negligible. 


(6) No change of phase takes place; all heat 
transferred is sensible heat. 


(7) There is equal heat transfer surface in each pass. 

The configuration is shown in Figure 3.1 where the five 
tube passes are designated with subscripts a, b, c, d and e. 
The temperature of the hot (shell side) fluid is indicated 
by upper case letters. For the cold (tube side) fluid, 
lower case letters are used. The subscript 1 always refers 
to the fluid inlet and the subscript 2 always refers to the 
fluid outlet. 

With Wy, and Coh designating mass flow (kg/sec) and 
specific heat (Joules/kg™°K) of hot fluid entering at T, and 


leaving at TI5 we define a Capacity tate ston sr iemierolae 


Ch = WhCon 
24 


Temperature 


| 
TT Y 
aul 
y 4 
| | 
| t 
| | V 
te 
| 
| 
{ | | 
de ree 
d | 
ted 
i | 
t | | 
“i ty 
tab 
| | 
t : | 
_dx | 
Xx 
0 a 


Figure 3.1 Three Counter Two Parallel Pass Configuration 
for Development of the Sought After 
Effectiveness Relationship 
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In similar fashion for the cold side (with We and Coc) 


entering at tj; and leaving at t9, we have 
Co = Welpe 


We then obtain an energy balance for the entire 


exchanger 
Cyn (11 = fo) = © oaaetre (3.5) 
Over the right hand side of the exchanger (Figure 3.1) 
Cyh(Ty - T) = Co(tg - te + tg = te + th - ta) Come 
and a differentiation gives 
ChdT = Ca(dte - dtqg + dt, - dtp + dtg) (Sip) 
Across dx, with a (m@/m), the surface per running meter 
of length of pass so that A = SaL is the total surface in 


the exchanger, we may write the heat transferred to the 


element dx in each cold pass. 


Codta = UVadxCl ata» (3.4a) 
Cadtp = -Uadx(T - tp) (3.4b) 
Cedt. = Uadx(T - te) (3.4c) 
Cedtq = “Uadx(T - tg) (3.4d) 
Cadte = Uadx(T = te) (3.4e) 
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Here it should be observed that due cognizance has been 
finem Ol the direction Of the flow in each cold fluid pass 
with respect to the positive sense of the length coordinate, 
x, and U is the overall heat transfer coefficient (W/m2—-°C). 


With eqs. (3.4) in eq. (3.3) 


or 
qT . = = oD 
ax — ny (oT mt a a th oe tq ae te) ( ° ) 
where 
Ua 
nN = 
h Ch 


is a sort of Ne, per unit length for the hot Side. 


Now differentiate eq. (3.5) 


q°T _ gat 7 dt, - dt, 7 dt. : dt, 7 dt. 
dx h “dx dx dx dx dx dx 


and with eqs. (3.4) substituted 


2 
a°T dT _ _ _ 
Bde TeTn’? ~ *a * *y - fe * fa” fe K8°9) 
where 
UVa 
an Oe 
c 


where again the resemblence of ne to Ney can be noted. 
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From eq. (3.2) we obtain 


- Tiger t, = atte eS Ue (3.7) 


1 2 b a Cc d 


and with eq. (3-7) put into eq. (3.6) 


2 

d“T dT _ _ h : a 

aye on), aaa noon a: a (T, 5h) to] 

dx e 
or 

ot _ g, (a ((R.- 1)T + T. - Rta] 

dx h dx coh C. Cc 1 c “2 

(3.8) 

where 

Re = Ce/Ch 


is the capacity rate ratio for the cold side. 


Notice that 








Con C. C. Cy ore C. 
and 
Hh. 7h ee Ce 
Cc Cc 


a capacity rate ratio for the hot side. Then, algebraic 


adjustment provides 
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—;——)T = m( 2 - ee) (33) 
Ry 1 / 


which is a linear, non-homogeneous, second order 
differential equation with constant coefficients having a 


complementary function 


5e (3.10) 


where C; and Co are arbitrary constants and where 


on 1 -R WZ 
= h ih Py h 
Sos) ess a ae a 5 [25n, oo ear 
_ on oy ae A ha = pe 
Ze 2 a 2 Ry 
h 
But 
RB eye a 
ieee la) oe h = 4 hy Bop we 
2 ae eo Cam h a) 
ny (C.) (Ua) R 
e @ 
so that 
= 9 
4,8, = °F + ah (Zeya Z i 2 ae 
) eae Ce 2 h Ry 
or 
Mh 
Ays% = 5 Care EA) Cone) 


Zo 


where 


h = (25 - 4R,(1 - Ree a (3.128 


Designate the particular integral as Tpi and by the 


method of undetermined coefficients let Tpj = P so that in 


eq- (9) 
1 - 1 
m ( a “ay p = m( a= - T,) 
h h 
This makes 
ie R 
7 = Zs h 
so that 
- ay 
pi 1 - R, (3. fae 
The general solution to eq. (3.9) is the sum of eqs. (3.10) 
and (3.2.13) 
Cia ax t, = RE 
T(x) = C,e : + Coe : oF oie a (3s 
h 


where the arbitary constants, C; and Co are evaluated from 
conditions at x = 0 and x = Ly At x = 90) T(x - 0) ] 153 
at x = L, T(x = L) = T,. When these are inserted, in turn, 
into eq. (3.14), one obtains a pair of linear algebraic 


equations in the unknowns C, and Co 


50) 


a. 
ul 
‘?) 


where Toi is given by eq. (3.13). 


It is only a matter of algebra to show Ena 


dob 
a (T, - Ti? - (T, - toa ve 
ii. .2. fue eo’ he 
al ay Ao 
e ——C 
and 
ah 
a (To ne (T, - om 
2 a,b dob 
e atl S. 


It is easy to see from eq. (3.1) that 


Ch (to = t,) 
R = = 
Be ON a 


so that 


to een) = eo) 


Use of this in eq. (3.13) shows that 


Bees oP nd 


pi en 


al: 
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G3. 154)) 


(3.15b) 


or 


by als 


Tot ae (3.16) 


indicating two alternative forms for Topi given by eqs. 
(3.13) 2nd (3 loo 
Insertion of eqs. (3.13) and (3.16) in eqs. (3. 15) ten 


Cy and C9 will yield after some algebra 


T, - to To = ty dob 
(= Sp yet 
h h- 
C, = : ——- (3.17a) 
1 2 
e -e 
and 
T. - ty a, Ty =_ to 
‘Tosmne are 
Cy = >t RP (3.17b) 
e -e 


Equation (3.14) is an expression for the hot side 
temperature at any location in the exchanger in terms of the 
extreme temperatures, t1, tg, T, and 19. 

Next take eq. (3.5) and set it equal to the derivative 


Of €q. (3524) noting sta te cic occ Tpi are all known 


constants. 

a a tC) — t. —*) ee 6 eae = C mis + C “25 

ge aS hb a ed be) eee aoCoe 
(3-1) 


At x = 0, where T = To, tg = t1, thy = te = the, tq = te = tde 


yy 


aT 


= ny (ST. = ty = 2trnan ZO =.) les a1Cy + AoCo G3219)) 


and if we subtract eq. (3.4a) from eq. (3.4c) we obtain 


C2 ome Ua oe 
on! Sealer: Teas = omy dx n dx 


which can be integrated using C3 as the constant of 


integration 


At x = 0 where t, = t, and te = tye 


ti - the = C3 


Or 
toe = t1 - Og 
In addition at x = L, tg = tap and te = teg so that 
“nb 
a eeced . 3° 
or 
- L 
ae tap ~ tea _ > : v3 
3 wah ab Ga 
S 


This gives a relationship between tgp and tne 


N 
C 


= 320 
ab Vag) ( ) 


B18) 


where No = noch can be considered as the total number of 
transfer units for the cold side. 
The same procedure may be employed to find tgqe- 


Subtract eq. (3.4c) from eq. (3.4e) 


Oe Sole 
e > 2a ae 
= ae = om dx n dx 
e 
and by integration 
=H x 
we = Ee = C4e 


At x = O where te the and te = tde 
tbe — tde = ©4 


or 


tde = the - C4 


Be 26 SIL els inles 1a teq and te = tg 


and again with noeL = No 


No 
CA = Cle deme = 


and 


N 
tae = the = acted saeco) cme (3.205 
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As before substract eq. (3.4b) from eq. (3.4d) 


at x = 0, th = type and tg = tge so that 


the — tde = C5 
and at x = L, tp = tap and tg = teg.- This shows that 
Neb 
Cab ~— tea = C5e 


and 
No 
teq = tab - (the - tdede (3.22) 
wmoere No = nolL. 
Return to eq. (3.18) and look at the conditions at x = L 
ere et, = th = tap, te = tq = teq, te = to and T = Tj- 


These conditions in eq. (3.18) give 


ny, (oT, - at. 
C3235) 


where again we remember that Cy; and Co are known constants. 
Now with four equations, eq. (3.20), eq. (3.21), 


eq. (3.22), and eq. (3.23) and four unknown variables we 
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can solve for the four variables, 


tab» the» 


ted, 


Cramer's rule will be used to solve for two of them; 


tde- 


and te.- 


toe 


First arrange eqs. (3.20) to (3.23) in standard form 


with the unknowns appearing in columns. 


+ “Ne 
tab “Hee 
2 “be 
tab ~ the 
tab 
where 8B = 1/2(5T, = ty 
ise 
-N 
ae. 
-1 
AX =B= 7 
e ¢ 
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The solution for the and 


; ' det Ane 
be 


det A 
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(3. 


(3. 


(3. 


and 


20) 


yi 


225 


23) 


where, Apa and Age are the matrices obtained by replacing 
the entries in the second and fourth columns of A by the 
entries in matrix B. 

Compute the determinant, det A of matrix A (the matrix 


of coefficients) 


=i 
ee es 0 
N 
Q =) ee. o4 
det A = . . 
1 e eC -1 = © 
i 0 1 0 
N -N 
det A = -(e © + 2e %) 


Since det A # 0, the system should have a unique solu- 


tion. Continue by computing det Ape and det Age- Thus 


—-N 
i t,e © of, 0 
N N 
0 toe & e 1 
det Ape = 
N 
1 O = ie 


ON N -N 
=weme (8 = Mee ia Ue ok Wye oS 


og 


and 


-N -N 
1. e . =I t,e : 
Ne Ne 
det Age = | 7 = v8 
N 
1 oe =i 0 
1 0 1 B 


2No . No 
= (e + 18 = Zt t,(e + 2e ) 


Therefore 
2N N -N 
dota e “(tae ot ene © + 2e5n@ 
Ucn ace Cae a a 
-(e + 2e ) 
or 
N 
€ “(2 =i 2to) 
ae = t = a (3.24) 
1 + 2e 
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det A (e + CC oe ee ce + 2e ) 
t _ de _ 74 1 
de det A No ms 
-(e + 2e ) 
or 
2Ne 
(e sue) i) = ato) 
Vae 7 Ln aan ee ee (3.25) 
e + 2e . 


ye! 


fimt@eme with eq. (3.24) and (3.25) in eq. (3.19) and some 


algebraic manipulation 


No -No 
Ge- 20,)(@e - er) 
a,C) + a,Cy = n, (5(T, - t,) + 2¢ ———+—_,, __—_ 
(ib 2e ~~) 
(3.26) 


Equation (3.26) is similar to the one derived by O'Hare 
[Ref. 1: eq. (22)]. At this point he continued to attempt a 
eoesed fOrm solution for €« as a function of R and Ny. It 
is to be noted, however, that he was unable to obtain the 
sOlution because of the presence of terms deriving from 
dimensional parameters such as ny, with dimensions m1 
[Ref. 1: pp. 46]. Since n, is present in eq. (3.26), it 
must be concluded that a further attempt at an analytic 
solution for the effectiveness of the 1-5 exchanger will 
suffer a similar fate and be also unachieveable. 

This conclusion should be obvious because, when dealing 
with an equation derived from n equations with ntl unknowns, 
One cannot create an additional equation by merely 
multiplying one of the n equations by a constant. This 
leads to a linearly dependent set which has no unique 
solution. 

Thus, after observing that O'Hare's attempt failed, this 
thesis, like O'Hare's, turns to a numerical analysis to 


achieve the requisite objective. 
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IV. NUMERICAL AND COMPUTER ANALYSIS 


As O'Hare discovered with the 1-3 heat exchanger 
[fRef. 1], an analytic solution for effectiveness as a 
function of the capacity rate ratio, R, and Ny, was also 
unattainable for the 1-5 exchanger. Hence, an alternate 
approach had to be taken and this involved utilizing a 
finite-difference thermal analyzer combined with a 
linearizing scheme for the solution of a set of nonlinear 
algebraic equations [Ref. 14]. This thermal analyzer has 
been employed to provide solutions for temperature that lead 
to effectiveness, e«, value for the 1-5:3C and 1-5:3P heat 


exchangers as a function of R and Ny,- 


A. THERMAL ANALYZER TVSSI 

The computer program TVSSI used in the development of 
the « = f(R, Ny) graphs for the 1=5 heat exchanger [Ketuauae 
Appendix A], was adapted from TVSS2 listed by Kern and Kraus 
[Ref. 14: Appendix C]. This adaptation consisted of modi- 
fying TVSS2 to compute solutions in the SI system and accept 
an input file specifically designed for tne 1-5 heat exchan- 
ger. This program also had to be adapted to the IBM 3033 AP 
system. The program utilizes the Cholesky decomposition 
scheme, and because of the linearization of the radiation 


terms (not used in this study) the program is iterative. 
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B. DEVELOPMENT OF THE 1-5:3C and 1-5:3P HEAT-EXCHANGER 
MODELS 


As indicated in the foregoing, TVSSI uses an input file 
to generate a solution for temperatures that will lead to 
the sought after effectiveness as a function of R and Ny. 
This input file is created by a program which requires an 
input of given capacity rates (hot and cold), overall 
coefficient, surface area and inlet stream temperatures (hot 
and cold) values. 

This required the development of a program which created 
such a file. O'Hare [Ref. lL: pp. 53-53] first developed a 
program using the 1-4 heat exchanger that would create an 
input file for TVSSI. He compared the effectiveness from 
the computer generated results and a known analytic solution 
of effectiveness developed by Kraus and Kern [Ref. 8]. Once 
O'Hare established confidence in his results, he proceeded 
to develop a program which would create an input file for 
the 1-3 exchanger. In this Study, a similar technique was 
used to develop a computer program for the 1-5:3C and 1-5:3P 
heat exchangers. These programs are called NTUS53C and 
NTU53P (see Appendices A and B) and, in them, the following 
parameters and techniques were utilized. 

1. 300 nodes were used. 


2. The initial temperature for the computer to begin the 
iterative process was set at 200°C. 


3. An eventual accuracy of .05 between the final and next 
to last iterations was used. 
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4. A radiation coefficient convergence factor of 0.66667 
between iterations was used. This is required by the 
program even if radiation is not considered as a 
node. 


o. The maximum number of iterations that the computer 
was allowed to perform was set at 12. 


6. A damping factor of .8 was set as an initial damping 
based on the number of non-linear terms in all of the 
node equations. 

With the values of Ch, Ce, U, A, Ti and ty, set, the 
input is set and an input file for TVSSI was generated. In 
this file all node equations and internode conductance 
values were determined. - This amieyeteeun specified the nodes 
that interact with each other and the methods by which the 
thermal interaction takes place such as in the forced 
convection and fluid flow modes. 

Both programs make use of the fact that each term 
in a node equation shows three things. The first is the 
node that is coupled for heat flow with the node in 
question. The second is the method of heat flow between the 
nodes. In this case,forced convection and fluid flow are 
used. Finally, the node equation shows the magnitude of the 
internode heat flow. Here all the pieces of information are 
collected and presented for use by TVSSI as an input file 
with all items in the proper format. 

It is the output values of both NTUS3C and NTUS3P (a 
typical result is shown in Appendix C) that are used by 


the thermal analyzer (as an input file) to determine 
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the temperatures Tg and tg for the specific set of given 


fierial Darameters C,, Co, R, A, T1, and tj. 


C. SCOPE OF COMPUTER ANALYSIS 

With the development of NTU53C and NTU53P computer 
programs, the files that they generated could then be used 
to solve for temperatures which in turn could yield effec- 
tiveness values for a given R and Ny,- Of course, it must 
be realized that many computer runs are required to generate 
enough data to ensure confidence in the results which cover 
a wide range of capacity rate ratios and Ny, values. 

To efficiently expedite the computer task, the Multiple 
Virtual System (MVS) with Job Entry Subsystem and Networking 
(JES3) was employed. The MVS coupled with JES3 is more 
commonly referred to as batch processing. Based on trial 
and error, it was determined that in order to build a solid 
data base, eleven different values for effectiveness were 
needed to best represent a particular value of R. Thus, a 
range of R from R = 0.01 to 1.0 was used. Due to the 
enormous amount of computer time required to solve a 300 
node problem, R increments of 0.05 were used for this study 
and this thesis reports on this basis. In all, 42 curves for 
both the 1-5:3C and 1-5:3P heat exchangers were required. 
With this change in the increments between R values 
(O'Hare's work showed that an increment in R of 0.01 was not 
necessary), a linear interpolation was found to be accepta- 
ble to less than 0.02 percent error. 
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TVSSI was slightly modified in accordance with the 
appropriate guidelines of the job control language (JCL) 
needed to run on the batch processing system. These 
modifications were few and were needed only at the beginning 
and end of TVSSI. The modified version of TVSSI has been 
called TVCOUNT with changes shown in [Ref. 1: Appendix F]. 
It was TVCOUNT that was then used to activate TVSSI. 

It also became necessary to modify the two input file 
programs NTUS3P and NTUS3C so that they needed to be 
compiled only once. They were nea loaded in a library 
file to be used when called by another program. New 
programs utilizing the batch system were written that could 
easily be loaded with the appropriate input data for a 
specific R value. These, which are referred to as "sister 
programs,'' were used to go from the library file to TVSSI 
and cause TVSSI to be executed eleven times under TVCOUNT 
covering the desired range of Ni, for 2 Specific A values 
The modified input files and their associated "sister 
execution programs are found in Appendices D through G. 

A flow chart of how all of the foregoing is accomplished 
is provided by O'Hare [Ref. 1: Figures 4.4 and 4.5]. It is 
noted that in these figures, TVSSI is referred to as TVSSIA 
through TVSSIV. These are the same programs as TVSSI but 
for bookkeeping purposes by the computer they are labeled A 


through V. 
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D. INTERPRETATION OF GRAPHS 

Once all the data had been collected and the 
effectiveness for a given R and Ny, values had been 
computed, plots of effectiveness as a function of Ny, for 
the entire range of R values were plotted for both parallel 
and counterflow heat exchangers. These plots are shown in 
Appendices H and I. 

It may be observed that for a particular value of R, the 
plot of effectiveness as a function of Ny, indicates that 
the 1:5-3C exchanger outperforms the 1-4, 1-3:2P and the 
1-5:3P exchangers but not the 1-3:2C exchanger. This may be 
explained by the quite apparent fact that the 1-5:3C 
exchanger has one more counterflow pass than the 1-4 
exchanger. Indeed, this argument can be extended to the 
fact that the 1-5:3C has one more parallel flow pass than 
ememi-3:2C. Figure 4.1. displays the graphical picture to 


Support this contention. 
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Figure 4.1 Comparison of Analytical (1-4) and Computer Results 
(1-3:2C and 1-3:2P) to Computer Results) (1—5:746 
and 1=5:3P) at Shao 
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V. POLYNOMIAL REGRESSION 


A. DEVELOPMENT OF POLYNOMIAL ‘EQUATIONS 

The empirical data obtained for the two 1-5 heat 
exchangers was designed to cover an extenSive range of R 
values varying from 0.01 to 1.0 in increments of 0.05. As 
discussed earlier in Section IV, the data obtained at a 
specific value of R is the computer evaluated result of the 
eifectiveness, for an associated N;, value. With this 
accomplished, it then becomes possible to graph separate 
eurves for each of the different R values as shown in 
Appendices H and I. Through a polynomial regression 
technique, as discussed in this section, it is also possible 
to develop implicit equations for the curves with € = fFCNeu; 
R). It is also apparent from an inspection of the graphical 
representation of the empirical data in Appendices H and I, 
that the curves conform to a high degree polynomial. 
However, further analytical investigation is needed to 
mecertain the exact order of the polynomial terms. This 
/Mavestigation will not only lead to the order of the 
polynomial, but to the specific equation for each curve. 

By use of polynomial regression, the least-squares 
method can be readily extended to best fit the data to the 


mth_gegree for the polynomial 


y = Ag + Ayx + Apx? + 226 Amxm (5.1) 
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with the error defined 


64 °> Ui a 


where Di represents eine 


to xj, x; being freevor 


The objective is to 


by 


empirical data value corresponding 
error. 


minimize the sum of the squares of 


the residuals, S,, 
m m Z 
» Za Z m 
Because at a minimum, the partial derivatives 3S,/9Ao, 
9Sy~/9A, .-- 9Sy/9A, vanish, after taking the derivative of 
Sy with respect to each of the coefficients of the pole 
nomial, it can be seen that 
aS 
iy me Z 7 : Zn z m 
aS 
‘ge = 2 2 
3A, = 0 = 2 x, (D; = Ag A,X; AgX; eee Ad ) 
aS 
ce oe 2 7 _ i 2 zi m 
BAS = 0 2d a8 (D, Ay A,X, AgX; eee A mo ) 
35, m 2 m 
TA = 0 = 2 1. Xs (D; - Ap - A,X; ~ AgX, — «ee -~ A x. ) 


Then by dividing by -2 and rearranging we obtain 
AgN + Aj = xq + A1xi- tees + Ay ye AU = 2 De 
Some + Ay © xy" + Ap © xy2 +... + Ap © x4Mtl = 5 x,D; 


Pex. = + Ay Fy xs + Ao £ x:* +... + AY 5 x-Mt2 = g x. 2p. 
0 7 aL i 2 a m i i °¥i 


eek, + Ay £ xiMtl + Ap xy xyMt2 + 1... + AYE xyoM = 5 x, Dj 


where all summations are from i=1 through n. All of the 
foregoing mt+l equations are linear and have mt+l unknowns: 
eee, 49, -.. A,- The coefficients of the unknowns can be 
calculated directly from the observed data. Thus, the pro- 
blem of determining a least-squares polynomial of degree m 


is equivalent to solving a system of m+tl simultaneous linear 


equations. Putting the equations in matrix form yields 
N oxy I ox4? Tee 2 ee xe y 
” Xq > x42 2, x43 yy ie ree xqmrl 2 
Mee’ 6 KAS Z x44 ees tea alee a= le): 
» et) } x mel y x4 mr2 uy x4 mts ee) Aes «42m 2» 





[Ref. 15: pp. 302-309 and Ref. 16: 468-474]. 
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From this point on, one finds that it is best to use a 
computer to assist in solving the simultaneous equations and 
this will also help alleviate any ill-conditioning that may 
otherwise occur. An existing "curvefit" program available 
through NON-IMSL [Ref. 16] was used although some modifica- 
tions were made to the original program to best accommodate 
the goals of this work. 

To determine the order of polynomial that should 
eventually be used, one increases the degree of the 
approximating polynomial as long as there is a statistically 
Significant decrease in the variance a2, which is computed 


by 


ON Naan etl (5-3) 


In otherwords, the selection of the optimum degree poly- 
nomial is contingent upon a decreaSing variance and once the 
variance begins to increase, the degree of the polynomial 
becomes too high. For all cases, it was found that the 
e - Ny, developed curves are of the oth order. 

As shown in Figures 5.1 and 5.2 the computed values of 
effectiveness vs. Ny, for R = 0.1, 0.5 and 1.0 for both 
flow arrangements, (1-5:3P) and (1-5:3C), have been graphed 
and fitted by a Sth degree polynomial. Because all computed 
values for effectiveness follow a predictable trend, 


only a sample of the data covering the whole range of 
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values of R have been shown. It is clear that the graphical 
interpretation strongly backs what is known numerically 
from the polynomial regression technique. Thus, the 


relationship for e« = f(Nt,, R) can be found explicitly from 
e = AsNy,° + AgNty* + AgNty? + AoNty? + AiNty + Ag (5-4) 


waere the coefficients As, A4g, A3, Ag, Aj, and Ag relating 
to a specific value of R are found in Tables 1 and 2 for the 
lme=o7oC) and (1-5:3P) configurations. 

In order to better utilize the coefficients found in 
Tables 1 and 2, an interactive program has been developed 
uSing equation (5.4) and the 1-5 heat exchanger polynomial 
approximation coefficients. This program was written for 
the IBM 3033 using the FORTRAN VS compiler (See Appendix J). 
By entering the. values for NTU, R and the type of heat 
exchanger the effectiveness values are readily obtained 
directly or by linear interpolation. 

An example of how to work with the effectiveness by 


using the Equation (5.4) follows. 


B. NUMERICAL EXAMPLE 
1. Statement of Problem 
A 1-5:3C heat exchanger uses sea water (Cp = 4.045 
kJ/kg-°C) at an initial temperature of 20°C to cool engine 
oil, flowing at 200 kg/s from 60°C to 40°C. The sea water 


flow rate is 80 kg/s. If the overall heat transfer 
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coefficient is 40 W/m2°C, how much surface area is reqiiieet 
in the exchanger? 
2. Solu Eroen 

At an average oil temperature of 50°C, the specific 


heat is Cy = 2.004 kJ/kg-°C and an energy balance gives 
CWCpAT) sea water = CWCpAT engine oil 
(80 kg/s)(4.045 kJ/kg-°C)(tg - 20)°C 
= (200 kg/s)(2.004 kJ/kg-°C)(60 - 40)°C 


= 90 + SE eee = 44.8°C 


2 e 


The number of transfer units is given by Nty = 
UA/Cmin: To determine Cyj, we must compare the products of 


mass-flow rate and specific heat for sea water and oil. 


(WCp)sw = (80 kg/s)(4.045 kJ/kg~°C) 


323.6 kW/°C 


(WCy)oil (200 kg/s)(2.004 kJ/kg-°C) 


400.8 kW/°C 


Therefore Cyjyn = 323.6 kW/°C. 
The effectiveness is a function of R and Nyy. The 


capacity rate ratiomenh, is 


o4 








Because C, = Cnin, the required effectiveness is 


to =. 


1 44.8 = 20 | 
=e tie OUe= 20 


We must now go to Table 1 for the 1-5:3C arrangement with 


R = 0.8 and R = 0.85 and find the coefficients 


R = 0.8 R = 0.85 
Ag = 3.63419 x 1079 Ag = 3.82659 x 1073 
Ay = 0.93822 A, = 0.93407 
Aj = -0.67053 Ao = -0.67841 
On 275 ai Ag = 0.28109 
Aqg = -0.60558 x 1071 Ag = -0.62088 x 1071 
As = 0.54393 x 1072 As = 0.55850 x 107? 


A trial and error solution (assuming values of Nyy and by 


interpolating) gives a result 


Ney = 2.352 


for which 


A = Ney Cmin/U 


19.03 m2 


ne 
il 


ey) 


3. Observations 

The primary observation made here is that by using 
the Sth order polynomial equation with the appropriate 
coefficients found in Table 1 or 2, an accurate value for 
effectiveness can be found thus allowing one to solve for 
5 Gide variety of heat exchanger parameters. 

By comparing the value for effectiveness just 
computed with that of the 1-5:3P, 1-3:2P [Ref. 1] and 1-4 to 
1-12 even pass exchangers [Ref. 8], the 1-5:3C has a higher 
effectiveness than all (see Table 3). 

From this observation it is evident that the 1-5:3C 
exchanger outperforms not only the 1-4 to 1-12 even pass 
arrangements but also its counterpart, the 1-5:3P exchanger 


for these given values by at least 0.614%. 
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COMPARISON OF EXAMPLE RESULTS 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The objective of this thesis was to provide effective- 
ness values for the 1-5:3C and 1-5:3P heat exchangers. 
Although these objectives were analytically unachievable, 
they became attainable through a fifth order polynomial 
approximation from numerical data. The data obtained covers 
a wide range of capacity rate mateters for values of Ny, from 
O to 3.25. With the knowledge of the 1-5 heat exchanger 
type (1-5:3C or 1-5:3P), size (surface area), the overall 
heat transfer coefficient and the fluid flow rates, the 
Ney and capacity rate ratios may be computed and the 
exchanger effectiveness may then be determined from the 
appropriate coefficient tables or the effectiveness program 
provided. 

It is indeed obvious that the effectiveness parameter is 
most useful in performance estimation and analysis for all 
types of heat exchangers. The work done here pertaining to 
the 1-5 exchanger is very helpful in constructing a 
comprehensive picture of the performance attainable for the 
1-2 consecutively through the 1-6 for a diverse range of 


known conditions. 


B. RECOMMENDATIONS 

With proof in this thesis, and that of O'Hare [Ref. 1], 
that an odd number of tube passes per shell has a slightly 
better effectiveness when more than half the tube fluid 
flows in counterflow to the shell fluid, it is recommended 
that a feasibility study be conducted with regard to the 
structural and therma 1 problems that may arise in 


manufacturing and designing such heat exchangers. 
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